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Abstract 
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SIMULATION OF FLIGHT MANEUVER-LOAD DISTRIBUTIONS BY UTILIZING 
STATIONARY, NON-GAUSSIAN RANDOM LOAD HISTORIES* 


By Herbert A. Leybold 
Langley Research Center 

SUMMARY 

The basic load history used to conduct random-load fatigue tests is, in general, 
Gaussian in nature and its peak distributions, both maximum and minimum, are inherently 
symmetrical about its mean. This symmetry is a characteristic of gust loading. How- 
ever, for maneuver loads, maximum and minimum peak distributions are asymmetric. 

Thus the basic Gaussian random load history does not afford an accurate simulation of 
the aircraft maneuver -load environment. The purpose of the present investigation was 
to determine the feasibility of utilizing some stationary, non-Gaussian random load his- 
tories to simulate this environment. 

Five stationary, non-Gaussian random load histories were generated with the aid of 
a digital computer. The peak distributions of the generated random histories were deter- 
mined and compared with maneuver -load flight data. The results indicated that non- 
Gaussian random load histories could be used to simulate slightly asymmetric aircraft 
maneuver-load spectra (i.e., spectra of transport aircraft), whereas highly asymmetric 
aircraft maneuver-load spectra (i.e., spectra of fighter aircraft) could not be simulated 
as well with these generated random histories. 

INTRODUCTION 

In recent years, random-load fatigue testing of aircraft structural parts was pro- 
posed. A major factor which influenced this approach was the realization that maneuver 
loads vary randomly with time and that, under the existing state of the art, fatigue life 
could not be predicted in an accurate, deterministic fashion. Currently, with more sophis- 
ticated fatigue testing machines available (machines capable of applying complex load his- 
tories which more nearly simulate actual flight experience), random -load fatigue tests 
have become rather common. 

*This paper is based in part upon a thesis entitled "Simulation of Service Load 
Distributions Utilizing Stationary, Non-Gaussian Random Time Histories” submitted 
in partial fulfillment of the requirements for the degree of Doctor of Philosophy in 
Mechanical Engineering, North Carolina State University, Raleigh, North Carolina, 1969. 



The random load history most widely used for random -load fatigue testing is 
Gaussian in nature. The Gaussian random load history is so called because its instan- 
taneous loads are distributed normally; this distribution of instantaneous loads can be 
characterized by the well-known Gaussian probability density function. But the peak 
loads, which are considered to be the significant fatigue -inducing parameter of a random 
load history, are not Gaussian in nature. Rice (ref. 1) has developed analytical methods 
for predicting the peak -load distributions of Gaussian random load histories. Both max- 
imum and minimum peak-load distributions are symmetrical about the mean of the random 
load history. This symmetry is a characteristic of gust loading. However, the maximum 
and minimum peak loads from maneuvers are not symmetrical about the mean. For 
maneuver loads, these distributions are asymmetric, and thus the basic Gaussian random 
load history does not simulate the maneuver -load flight environment. The purpose of the 
present investigation is to determine if certain stationary, non-Gaussian random load his- 
tories can simulate the aircraft maneuver-load flight environment. 

Since no analytical methods are available for determining the maximum and mini- 
mum peak distributions of non-Gaussian random load histories, these distributions were 
determined numerically. A random -number generator was used to generate random num- 
bers having a constant density function and having limits of 0 and 1. These random num- 
bers were then transformed so that the probability density function of a discrete random 
load history composed of these random numbers had one of the following non-Gaussian 
distributions: Poisson, binomial, log-normal, Weibull, and exponential. These five forms 
of the density function were arbitrarily selected because they are widely discussed in the 
literature and are mathematically describable. The transformed random histories were 
analyzed to determine their peak statistics and were compared with peak statistics encoun- 
tered in flight. 


SYMBOLS 


E( ) expected value 

F(x) cumulative probability density function 

f(x) probability density function 

g acceleration due to gravity, 9.8 m/sec^ 

n number of events (positive integer) 

p probability of success 
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q 


probability of failure 


RMS 

At 

x 

<y> 

a,/3 

r( ) 

e 


M 

M log 

a 


a log 


root mean square 
time increment 
random variable 

set of uniformly distributed random numbers 
parameters of Weibull probability density function 
gamma function 

parameter of exponential probability density function 

parameter of Poisson probability density function 

mean of logarithms 

standard deviation 

standard deviation of logarithms 

variance 

GENERATION AND ANALYSIS OF RANDOM LOAD HISTORIES 


Uniformly distributed random numbers were generated digitally and then trans- 
formed to have a nonuniform distribution. The set of uniformly distributed random num- 
bers was transformed into sets having Poisson, binomial, log-normal, Weibull, and expo- 
nential density functions. These functions were chosen because they are common, their 
characteristics are well known, and they are analytic. The transformed numbers are 
considered to be the instantaneous values of load at discrete time intervals and therefore 
form a load history (see fig. 1). The final step was to study the peaks in the five load 
histories generated. These peaks are not distributed according to the distribution of the 
set of transformed random numbers; for example, a load history with instantaneous val- 
ues that are Poisson distributed has maximum and minimum peaks that are not Poisson 
distributed. Figure 1 shows the differences between what is meant by a peak (maximum 
or minimum) and the instantaneous value of the load. 
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A set of random numbers {y} bounded by the limits 0 and 1 and having a constant 
probability density function were generated with a standard random -number generator. 
This set of uniformly distributed numbers was then transformed into five new sets; each 
new set was distributed according to one of the five previously mentioned non-Gaussian 
probability density functions. The numbers, when drawn from a particular set, represent 
the instantaneous values of the load history; obviously, these instantaneous values of the 
history are distributed the same as the numbers of the set. 

The set of uniformly distributed random numbers was transformed by 

F(x) = C* f(y)dy (1) 

oo 

where f(y) is the specified density function. The cumulative probability density func- 
tion F(x) is bounded by the limits 0 and 1, where x can take on values of -«° to +°°. 
The random numbers from the set {y} having a constant probability density distribution 
and having limits of 0 and 1 were substituted for F(x) in equation (1). Equation (1) could 
then be evaluated for x (the transformed random number). For this investigation, each 
random history generated contained 10 000 random numbers. 

Two of the probability density functions considered (Poisson and binomial) were 
discrete functions, and hence the integral in equation (1) was changed to a summation to 
determine F(x); that is, 


A. 

F(x) = £f(y) 


( 2 ) 


For the Poisson distribution, the following recurrence formula was used to compute F(x): 


f(x + 1) = - z f (x) (x = 0, 1, 2, . . ., °°; M > 0) (3) 

X T X 

Similarly, for the binomial distribution, the recurrence formula used was 

f(x + 1) = f(x) (x = 0, 1, 2, . . ., n; 0 < p < 1; 0 < q < 1) (4) 


The density function of the log-normal distribution was 

\ 2 ~] 


f(x) = - 


-log^ 


exp 


-(logg x - Miog) 


2(7 


log 


(x > 0; <T log > 0; -oo < M log < °°) (5) 


where x was evaluated with a numerical approximation of the integral in equation (1). 
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For the other two density functions, x was determined by directly integrating 
equation (1). For the Weibull distribution, 


f(x) = a/3x a ~lex p(-0x a ) 


and 


F(x) = 1 - exp(-/3x a ) 
Solving for x in equation (7) yielded 



(x > 0; /3 > 0; a > 0) (6) 

( 7 ) 


( 8 ) 


For the exponential distribution, 


f(x) = e exp(-#x) 


(9 > 0; x > 0) (9) 


and 


F(x) = 1 - exp(-#x) (10) 

Solving for x in equation (10) gave 

x = “ \ lQ ge & " f ( x 3 (11) 

It should be noted that the exponential distribution is a special case of the Weibull distri- 
bution when a = 1. 

The first 100 transformed random numbers for each of the five different load his- 
tories generated are plotted in sequence at equal intervals of time in figures 2 to 6. The 
mean, the root mean square, and the standard deviation are shown for each history. 

The mean indicates the average or expected value of the random variable x and 
is mathematically expressed as 



xf(x)dx 


( 12 ) 


The root mean square is a measure of the amplitude of a time history and is the square 
root of the mean of the squares of a random variable x. The mean square (that is, the 
mean of the squares) is mathematically expressed as 
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( 13 ) 


e(x 2) = f x2f(x)dx 

oo 

The standard deviation is a measure of the dispersion from the mean and is the square 
root of the variance of a random variable x. The variance is defined mathematically as 

ct2 = Ejjx - E(xj] 2 J (14) 

which reduces to 

2 

ct 2 = E(x2) - [ E (x)] (15) 

All three properties were calculated for each of the five random load histories gen- 
erated. The results of these calculations are shown in table 1. These same properties 
for the selected functions were calculated on the basis of the analytic equations summa- 
rized in table 2. The resulting theoretical values showed good agreement with the values 
calculated for the generated histories, as can be seen in table 1. 

COMPARISON OF GENERATED AND THEORETICAL HISTOGRAMS 
OF INSTANTANEOUS VALUES 

The frequency of occurrence of the instantaneous loads was determined (see fig. 7(a)) 
for each of the five random load histories generated, and these frequencies are plotted as 
histograms in figures 2 to 6. The distributions of the generated numbers are compared 
with theoretical histograms. The theoretical histograms were obtained by first taking an 
incremental area under the particular probability -density -function curve and then multi- 
plying the incremental area by 10 000 so that the frequency of occurrence has the same 
base as the generated numbers. Figures 2 to 6 indicate that excellent agreement was 
obtained. In addition, a check on the randomness of the number generator was performed 
with a chi-square test. In all cases, the chi-square value is inside the accepted region 
specified by the 95-percent confidence level. 

Some of the parameters selected for the binomial, Weibull, and Poisson density 
functions tended to make them look Gaussian. This result was to be expected because the 
normal or Gaussian density function is frequently used to approximate these other func- 
tions. For example, the normal function approximates the binomial function quite well 
when the binomial parameters np and n(l - p) are both greater than 5, as verified in 
figure 5(b). 
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DETERMINATION OF PEAK AND CUMULATIVE PEAK DISTRIBUTIONS 
FOR BOTH MAXIMUM AND MINIMUM PEAK VALUES 

Peak values, both maximum and minimum, were counted in equal intervals as shown 
in figure 7(b). The frequency of occurrence of these peak values, for each of the five 
random histories generated, are plotted on the left-hand side of figures 8 to 12. The 
maximum and minimum peak distributions intersect at the mean of the random history. 
Moreover, for Gaussian random histories, both maximum and minimum peak distributions 
are approximately normal in shape and are mirror images of each other. As mentioned 
previously, some distributions of instantaneous values show a trend toward being normal. 
Therefore, it could be anticipated that their maximum and minimum peak values are 
approximately normally distributed and symmetrical about the mean (see, for example, 
figs. 5(b) and 11(b)). 

Cumulative peak distributions were also determined and are plotted on the right- 
hand side of figures 8 to 12. Some of these distributions are subsequently compared with 
actual maneuver -load data from aircraft in flight. 

COMPARISON OF CUMULATIVE PEAK DISTRIBUTIONS 
WITH MANEUVER-LOAD DATA 

Actual maneuver -load data are generally presented on the basis of either the incre- 
mental acceleration (in g units) or the percent limit load. The data are usually plotted 
on semi -log paper with the cumulative number of exceedances plotted on the log scale and 
the incremental acceleration or percent limit load on the linear scale. The digitally gen- 
erated data discussed previously have been converted to these scales for purposes of 
comparison with actual maneuver -load data. In order to cover a wide range of asymmetry 
of actual loading spectra, two different types of spectra, one for a fighter aircraft and one 
for a transport aircraft, were selected for comparison. The fighter spectrum is highly 
asymmetric whereas the transport spectrum is only slightly asymmetric. For the range 
of parameters covered in this investigation, all peak distributions generated were com- 
pared with these two types of spectra. 

Figure 13 presents the operational maneuver-load spectrum for transport -type air- 
craft as obtained from reference 2. This spectrum is compared with the digitally gener- 
ated data which have a log-normal probability density function. Good agreement between 
the digitally simulated and actual maneuver -load data was obtained. The simulated data 
resulted in more occurrences than did the actual flight data, particularly at the lower 
levels of incremental acceleration. 
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Figure 14 presents the operational maneuver -load spectrum for fighter -type air- 
craft as obtained from reference 3. This spectrum is compared with the digitally gener- 
ated data which have a Weibull probability density function. Fair agreement between the 
digitally simulated and actual maneuver -load data was obtained, the best agreement being 
achieved for the positive half of the spectrum. 

CONCLUSIONS 

Five stationary, non-Gaussian random histories have been generated with the aid of 
a digital computer. Distributions of maximum and minimum peaks as well as probability 
density functions have been obtained for each load history. The following conclusions 
have been drawn from an analysis of the generated digital data: 

(1) When the probability density function of a random history approached a Gaussian 
distribution, the maximum and minimum peak distributions of that random history were 
mirror images of each other about its mean and were distributed essentially normally. 

(2) When the probability density function of a random history was decidedly non- 
Gaussian in shape, the maximum and minimum peak distributions of that random history 
were asymmetric about its mean and were also non-Gaussian in shape. 

(3) Of the five random histories generated, the one having a log-normal distribution 
compared most favorably with the transport -type spectrum whereas the one having a 
Weibull distribution compared most favorably with the fighter -type spectrum. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 22, 1971. 
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TABLE 1.- PROPERTIES OF GENERATED AND THEORETICAL DENSITY FUNCTIONS 


Density 

Parameters 

Mean, E(x) 

Root mean square, RMS 

Standard deviation, a 

function 

Generated 

Theoretical 

Generated 

Theoretical 

Generated 

Theoretical 

Weibull 

a = 2; /3 = 2 

0.6207 

0.6267 

0.7004 

0.7071 

0.3245 

0.3276 


a = 5; /3 = 5 

.6631 

.6655 

.6801 

.6827 

.1511 

.1524 

Exponential 

0=1 

0.9811 

1.0 

1.3917 

1.4142 

0.9870 

1.0 


0 = 2 

.4906 

.5 

.6959 

.7071 

.4935 

.5 

Poisson 

jx = 10 

9.948 

10.0 

10.436 

10.488 

3.130 

3.162 


p = 66 

65.867 

66.0 

66.365 

66.498 

8.032 

8.124 

Binomial 

n = 10; p = 0.1 

0.979 


1.357 

1.378 

0.941 

0.949 


n = 500; p = 0. 1 

49.892 



50.448 

6.626 

6.708 

Log-normal 

%g ‘ 0; '’log = <U 

1.0453 

1.0513 

1.0983 

1.1052 

0.3372 

0.3409 


•’log = °> T fog “ °- 05 

1.0214 

1.0253 

1.0468 

1.0513 

.2296 

.2322 


to 













TABLE 2.- SUMMARY OF PROPERTIES OF DENSITY FUNCTIONS 


Density 

function 

Mean, E(x) 

Mean square, e(x^) 

Variance, a 2 

Weibull 

, \l/a , N 

(?) r h') 

1 

(l\ r /a£ + 2 \ 

W V “ / 

/ \ 2/a 

V/ 3 / 1 

r fl 


• 

Exponential 4 

! 1 v 

i 

‘ 

2 

0 2 

1 

e 2 

1 *i 

i Poisson ju /x(p. +1) jx 

i j 

; Binomial : np 

j 

n(n - l)p2 + np 

np(l - p) 

Log-normal ! exp (d log + \ o? og ) 

exp (^Miog + ^ a iog) 

exp(2Mi og + o2 og )exp(c2 og )-l 





Figure 1.- Discrete random history showing maximum and minimum peaks. 


Relative load Relative load 


2.5 


° Generated 
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(a) a - 2; /3 = 2. 



(b) a - 5; /3 = 5. 

Figure 2.- Sequences of first 100 random numbers and histograms of first 10 000 random 
numbers for load histories having a Weibull probability density function. 
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Relative load Relative load 


10 


o Generated 
— Theoretical 



jj me — Frequency of occurrence 

(a) 0=1. 



Time — Frequency of occurrence 

(b) 0=2. 

Figure 3.- Sequences of first 100 random numbers and histograms of first 10 000 random 
numbers for load histories having an exponential probability density function. 
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Relative load Relative load 


25 r 



(a) n = 10. 



(b) m = 66. 

Figure 4.- Sequences of first 100 random numbers and histograms of first 10 000 random 
numbers for load histories having a Poisson probability density function. 
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Relative load Relative load 



Time — - Frequency of occurrence 

(a) n = 10; p = 0.1. 



Time — Frequency of occurrence 


(b) n = 500; p = 0. 1. 

Figure 5.- Sequences of first 100 random numbers and histograms of first 10 000 random 
numbers for load histories having a binomial probability density function. 
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Relative load Relative load 




Figure 6.- Sequences of first 100 random numbers and histograms of first 10 000 random 
numbers for load histories having a log-normal probability density function. 
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Frequency of 
occurrence in 
each interval 



0 

1 

0 
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3 

3 

4 

4 

0 
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(a) Frequency of occurrence of instantaneous loads. 
Figure 7.- Frequency of occurrence in prescribed intervals. 
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(b) Frequency of occurrence of maximum and m in i m um peaks. 
Figure 7.- Concluded. 
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(a) a = 2; /3 = 2. 
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o I0 3 
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A .8 

Relative load 


(b) a = 5; p = 5. 

Figure 8. - Peak and cumulative peak distributions for random histories having a 

Weibull probability density function. 
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(b) 0=2. 
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Figure 9.- Peak and cumulative peak distributions for random histories having an 

exponential probability density function. 
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(a) 11 = 10. 



Relative load 
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(b) ii = 66. 

Figure 10. - Peak and cumulative peak distributions for random histories having a 

Poisson probability density function. 
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(a) n = 10; p = 0.1. 
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(b) n = 500; p = 0.1. 

Figure 11.- Peak and cumulative peak distributions for random histories having a 

binomial probability density function. 
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Figure 12.- Peak and cumulative peak distributions for random histories having a 

log-normal probability density function. 
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Figure 13.- Log-normal digitally generated data compared with operational maneuver -load 
spectrum (for transport aircraft) from reference 2. 
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Percent of limit load 

Figure 14.- Weibull digitally generated data compared with operational maneuver -load 
spectrum (for fighter aircraft) from reference 3. 
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